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Abstract 

We propose a model for a spatially modulated collective state of supercon- 
ducting cuprates in which the electronic properties vary locally in space. In 
this model the regions of higher hole density (called stripes) are described as 
Luttinger liquids and the regions of lower density (antiferromagnetic ladders) 
as an interacting bosonic gas of d x 2_ y 2 hole pairs. We show that the transi- 
tion to the superconducting state is topological and driven by decay processes 
among these elementary excitations in the presence of vibrations. 
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The experimental evidence for charge and spin inhomogeneities in superconducting 
cuprates has been accumulating over the last few years!. It is believed that these inho- 
mogeneities are the result of strong competing electron-electron and electron-lattice inter- 
actions. The strong interactions can be hinted from the fact that cuprates are directly 
related to Mott insulators and not to metals. There is an emerging point of view that in- 
homogeneities may be at the heart of the superconducting problem and that an important 
part of the physics occurs in real space. This should be contrasted with a Fermi liquid 
description which is dominated by the Fermi surface. In this work we propose a model for 
superconductivity which requires a dual description in terms of real and momentum space. 

We assume the existence of two main length scales in the problem. At the short length 
scales the physics is dominated by the strong Coulomb forces associated with the Mott 
insulator!. These forces lead to the quantum confinement of holes!. An example of such 
confinement happens when two holes are injected into an antiferromagnet and there is a gain 
in magnetic and kinetic energy if they move together instead of separately. This leads to the 
formation of bound states of two holes which due to the same strong Coulomb forces and 
lattice symmetry have d x 2_ y 2 symmetry!. On the other hand, in the presence of magnetic 
defects such as anti-phase domain walls one-dimensional (ID) structures such as "stripes" 
can be stabilized due to the gain of single hole kinetic energyR Stripes are regions of larger 
carrier density which have properties analogous to the ones found in low-dimensional inter- 
acting Fermi systems such as Luttinger liquids!. In this case a long length scale associated 
with the inter-stripe distance emerges leading to a phase with long-range orientational order 
in the absence of translational order, that is, a state similar to a classical liquid crystal!. 
In fact, because the charge modulated state (CMS) is composed of electrons we expect this 
modulated phase to be quantum in nature, that is, a quantum liquid crystal. 

Let us consider a static CMS induced along the y direction in the 2D system with peri- 
odicity Na where a is the lattice spacing. We assume that the CMS is a superposition of 
a dilute homogeneous bosonic background of d x 2_ y 2 hole pairs in 2D and a set of ID Lut- 
tinger liquids. This picture is supported by density matrix renormalization group (DMRG) 



calculations of the t — J mo In fact, weak disorderly, spin-orbit effects and lattice 

distortions^ might stabilize density modulations even further. The CMS is such that bosons 
are continuously created and destroyed from the Luttinger liquids. In the absence of cou- 
pling between such excitations the system should be insulating even in the presence of the 
smallest amount of disorder: a Bose glass of pairs! and a charge density wave state (CDW)i 
on the stripes (we assume repulsive interactions only). Deformations of the amplitude of the 
density of the CMS are always energetically costly. Therefore, phase fluctuations are the low 
energy excitations in such systems. Moreover, charge neutrality implies that a CMS has to 
be coupled to the lattice. Therefore fluctuations of the CMS should appear directly into the 
phonon spectrum. We believe that this has been observed experimentally in essentially all 
the cuprateJH. The mechanism for superconductivity discussed in this work requires the 
coupling of stripe fermions via the exchange of bosons. One can think of this mechanism 
as the exchange of stripe "pieces". The simplest exchange mechanism is due to the decay 
of these bosons into fermionic degrees of freedom at the stripes. We find, however, that 
the decay process requires the production of lattice vibrations. As a natural consequence, 
superfluid vortices are coupled to dislocations of the CMS. Because the superfluid density is 
low, the interactions between topological defects ultimately determine the phase diagramS. 

We assume the existence of bosons with energy and local repulsion U described by 
creation (annihilation) operators (PjJ and fermions with energy dispersion ek and spin 
a (| or |) described by creation (annihilation) operators ip\ ia (V'k.o-)- When iV ^> 1 the 
direct exchange of fermions between stripes is suppressed and we can describe the fermions 
as isolated ID Luttinger liquids. In this limit it is convenient to label the fermion operators 
as ip n ,a,cr{x) where n is the stripe label, a = R,L refers to right and left moving fermions!. 
On the one hand, bound states of two holes are not eigenstates of the Luttinger liquid. On 
the other hand, collective sound waves of the Luttinger liquid are not elementary excitations 
of the weakly doped Mott insulator. Thus, the simplest process which couples these two 
types of excitations is a decay process. The decay produces a "shake" of the collective state 
emitting phonons. Conservation of linear and angular momentum requires the coupling 



between fermions and bosons to have the form (we use units such that % = ks = 1 and 
/3=1/T): 

H c = J2 ^(^ApPk^tV'-k+p-q,;^ + h.c. . (1) 

k,p,q 

where b p (tf p ) is the annihilation (creation) operators for vibrations of the collective state 
(including the lattice). Here A p — > when p — > is the stripe- vibration coupling constant 
(the momentum dependence is a consequence of the local coupling of the bosons with the 
stripes), Vk(q) is the boson-stripe coupling which, by conservation of angular momentum, 
has the d x 2_ y 2 symmetry. Notice that to destroy a boson is equivalent to create two electrons 
in the antiferromagnetic ladder. In fact, ([I]) is very similar to fermion-boson models which 
have appeared in the literature of superconductivity over the years@. It is straightforward to 
see from (p]) that the condensation of the bosons with simultaneous condensation of fermions 
into Cooper pairs with zero center of mass momentum requires superfluidity and lattice dis- 
tortions to occur at the same wave-vector. In particular, the largest coupling occurs at the 
edge of the Brillouin zone with wave- vector K = (0, ±ir/a) and thus perpendicular to the 
orientation of the charge modulation (that is, perpendicular to the "stripes"). In the case 
of the formation of domains with stripes running in orthogonal directions we also expect 
response at K* = (±7r/a, 0). If we forget about the phase fluctuations we would find a 
superfluid/superconducting phase with simultaneous dimerization of the CMS perpendicu- 
lar to its orientation. This result would be in agreement with the simultaneous appearance 
of the pseudogap in angle resolved photoemission experiments!^! (ARPES) and the phonon 
anomalies observed in neutron scattering^ exactly at K. We stress that although vibra- 
tions are not the mechanism of pairing (which is essentially produced the electron-electron 
interactions), they are essential for condensation. 

Because we are discussing a 2D system with low superfluid density the problem of su- 
perconductivity is decided by phase fluctuations. Consider the partition function of the 
problem which can be written in terms of an Euclidean Lagrangian L. We can separate 
the problem into fast and slow moving fields by linearizing the fermion dispersion close to 
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the stripe Fermi surface and studying the bosons and lattice distortions close to K. We 
disregard all the fast oscillating terms in the Euclidean action and rewrite the Lagrangian 
as L = L [P, b] + J2 n LsbPn] + L c where 



L [P,b] = \j d 2 r\p{v) (d T + E 



V 2 \ U 

* P( r ) + _AT 2 (r) 
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where E (cu ) is the energy of the bosons (vibration) at K, M B is the effective boson mass, 
My is the mass associated with the vibrations and N(r) is the local number of bosons, Ls 
is the Lagrangian for interacting fermions in mi and 



L c = \Y^ I dxV(x,n) P{x,n)b{x,n) + P(x,n)b{x,n)] [^ n ,R^(x)^ n ,L,i( x ) 

^ n J 

+ l/} n ,L,l(x)lp n ,R^{x)] + C.C.} 



(3) 

is the boson-stripe-vibration coupling. We now redefine fields as P = y/d^e~ lipB and b = 
■^f<jye +lipv and use the bosonization technique! in order to get L = J2n LiD[9 n ,p, 4>n,s] + 
J2 a =B,v L a [cr a ,<p a ] + L [a a ,(pa,9 P ,<fis] where L lD [9 n ^,4> nyS ] is the action for ID bosons (9 p 
and <p s are the charge and spin bosonic modes at the stripes)§ and 



d 2 r 
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ia a d T ip a + E a a a + ^-a 2 a + ^-(VcrJ 2 



(To. 



2M< 



"(V^a) 



(4) 



with Eb = E , Ey = ujq, Ub = U and Uy = 0. Finally, 

L c = V / ^K(x,nAra) v ^i^e- i( ^ (x ' nWa)+ ^ (x ' nWa) -^"' p{a;)) cos(v / 270„, s (x)) + c.c. . (5) 
~ J 2~na 

We can gauge away the phase fields and define a new bosonic field 9 C such that 9 njC (x) = 
9 n , P (x) — (ip B (x,nNa) + (fy(x, nNa))/\TM so that L c simplifies to 

<h ' _ ~ (6) 



/qx 
— V(x, nNa)y/a B <Jv cos(V2ir9 n c (x)) cos(V2ir<f) ntS (x)) . 
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Because of the gauge transformation a new term appears in the Lagrangian: 



L /M = S / dx 



(d T (ip B (x, n) + (p v (x, n))f + v 2 {d x ((p B (x, n) + tp v (x, n))f 



+ 



2itv r 



d T {ip B {x, n) + ip v (x, n))d T 9 n)C (x) + v 2 d x (ip B (x, n) + (p v {x, n))d x 9 n;C (x) \ . (7) 



Here K p is the Luttinger parameter and v p the charge velocity. We can study the phase 
fluctuations around the saddle point where: a a = cr°, (p a — with a = B, V. This just 
gives the mean field solution which will be studied elsewhere^. This mean field solution 
produces a characteristic temperature T* below which superfluidity and lattice distortions 
start to be simultaneously generated at K. This temperature, however, is not the actual 
transition temperature since the phases of the order parameters fluctuate strongly. T* 
is a crossover temperature which we associate with the so-called pseudogap phase of the 
cuprate Around the saddle point we can easily show that Lq in (f|) is a relevant operator 
in a renormalization group sense even in the presence of strong repulsion. Thus, even if 
isolated stripes are deep inside of a CDW state the coupling V drives the system to a 
superconducting state. Amplitude fluctuations around the saddle point are therefore massive 
and produce a spin gap in the system. Integrating over amplitude fluctuations up to second 
order leads to an Euclidean action of the form: 

s = 77 / dr [ d2r \ J2 i°ad T (p a (r, t) + ^ (d T cp a (r, r)f + ^j- (V<^ Q (r, r)f 
J2 <9 r ^ Q (r,r) +v 2 l J2 d x (p a (r,r) 



+ K > 



2iTVpNa 



,a=B,V / \a=B,V 



(8) 



where the factor of 1/N appears due to the coarse-graining of the fields in the direction 
perpendicular to the stripes. Here a B is the total density of electrons in the system, Kb is 
the charge compressibility (the superfluid velocity is given by cb = 4:(Xb/(kbMb)), ~oy = 
and is the lattice compressibility (the sound velocity is cy = 4<Ty /(k v M v )). The key 
point about (§) is that the topological excitations of the superfluid state (the vortices) are 
directly coupled to the topological excitations of the dimerized CMS (the dislocations loops) 
via the stripes. 

Let us consider now the problem at < T < T* so that the amplitude of the order 
parameters are well developed but phase coherence has not been established yet. At finite 
temperatures we can disregard the time derivatives in (|S]) and rewrite 
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Z J \ a=B,V 1V±a Z/UVU L =B y 



(9) 



which describes two 2D XY models coupled along the x direction (vf is the stripe Fermi 
velocity). In what follows we will study the problem of phase coherence due to the renor- 
malization of one of the phase modes by the gaussian modes of the other. This procedure 
works when the coupling is very weak, that is, N 1. Because of the quadratic nature of 
the action it is a simple exercise to show that the problem is described by a 2D XY model 
with renormalized stiffness. In the superffuid case the stiffness is given by 



Pb = (Eb + E\ 



1 E v j E B 



(10) 



'l + E c (E B l + Ey 1 ] 

where E a = a a /M a and Ec = vf/ '(2ttNo). An analogous result is obtained for the rigidity of 
the vibrations with the replacement of of Eb by Ey and vice-versa. Therefore, we conclude 
that the transition from the ordered to the disordered phase of the superfluid or/and the 
dimerized CMS is due to the unbinding of the topological excitations: vortex-antivortex 
pairs in the case of the superfluid and dislocation loops of the dimerized CMS. In fact these 
topological excitations are coupled since a dislocation of the dimerized CMS can potentially 
liberate vortex-antivortex pairs. The transition temperature to the ordered phase can now 
be estimated directly from (|10"D as 

Tc,a ~ \pa (11) 

in accord with the Uemura relation^. An interesting consequence of our calculations is 
that if the phonon at K is optical (dispersionless) then Ey — > and py — > 0. In this 
case the transition to the CMS is driven to zero temperature! Thus, while the system 
becomes a superfluid-superconductor at T^s, at any finite temperature dislocations loops of 
the dimerized CMS will render T c y = 0. 

In summary, we have proposed a 2D model for superconducting order in cuprates 
which involve a spatially modulated collective state which contains d x 2_ y 2 bosons and stripe 



fermions. We show that the coupling between these particles leads to superconductivity and 
simultaneous condensation of bosons and phonons at finite wave-vector. Although conden- 
sation starts at high temperature T*, phase coherence is only attained at low temperatures 
due to proliferation of vortex-antivortex pairs of the superfluid and dislocation loops of a 
dimerized CMS. 
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